UNSTEADY  CONJUGATE  HEAT  EXCHANGE  OF 
OIL-GAS  PIPE  LINE  IN  SOIL 

A.M.  Pavljuchenko  and  A.N.  Popkov 


Institute  of  Theoretical  and  Applied  Mechanics  SB  RAS, 
630090  Novosibirsk,  Russia 

Introduction 


The  problem  of  conjugate  heat  exchange  of  main  gas-oil  pipe  lines  is  of  special 
importance  for  pipelines  laying  in  the  permafrost  regions,  because  the  energy  saving,  reliability 
and  ecology  problems  ought  to  be  solved.  In  this  work,  this  complex  mathematical  task  without 
allowance  for  a  movable  thawing  bound  was  numerically  simulated. 

The  calculating  results  of  conjugate  unstable  heat  exchange  of  a  “hot”  product-pipeline 
with  a  uniform  soil  were  presented  at  the  following  assumptions: 

-  The  flow  regime  in  a  tube  is  turbulent  and  unsteady,  the  liquid  is  incompressed,  and 
the  physical  properties  are  constant. 

-  Heat  exchange  along  the  pipeline  axis  is  negligibly  small  in  comparison  with  the  heat 
exchange  normal  to  the  pipeline  axis. 

-  The  soil  in  which  the  product-pipeline  is  embedded  is  a  uniform  massif  with  known 
thermal-physical  properties. 

-  The  pipe  wall  and  its  insulation  is  negligible  thin,  and  their  thermal  resistance  is 
modeled  by  a  certain  coefficient  in  the  boundary  condition  at  the  pipe  and  soil 
junction. 

The  task  for  the  pipe  at  a  flow  hydraulic  stabilizing  comes  to  solution  of  the  energy 
unstationary  equation: 
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where  T  is  a  liquid  temperature;  r  is  a  time;  a  ,  Uj  are  the  coefficients  of  a  temperature 


conductivity  of  the  laminar  and  turbulent  flow  regimes  accordingly;  IE  is  a  velocity  axis 
component; 

Equation  (1)  describes  the  unsteady  heat  exchange  of  liquid  and  the  pipeline  wall  at  the 
turbulent  flow  regime,  and  the  velocity  field  is  described  by  the  Reichard  dependence  [1]  and 
[2].  Equation  (1)  is  a  hyperbolic-parabolic  (ultra -parabolic)  that  are  considered  at  difference 
scheme  construction. 

The  temperature  distribution  in  soil  at  undersoil  pipe  laying  taking  into  account  a  day 
surface  is  described  by  a  non-stationary,  two-dimensional  equation  of  the  thermal  conductivity. 
In  such  case,  the  longitudinal  heat  flows  are  neglected  in  comparison  with  the  transversal  ones. 
Thus,  a  three-dimensional  problem  of  the  thermal  conductivity  is  split  in  a  great  number  of  two- 
dimensional  problems  connected  by  boundary  conditions  on  a  pipe  wall: 
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t  is  a  soil  temperature;  is  the  coefficient  of  the  soil  thermal  conductivity;  x  ,  y  are  the  co¬ 
ordinates  of  the  cross  section;  x  is  directed  along  the  soil  bound  and  the  day  surface  from  a 
symmetry  plane;  axis  y  -  lies  in  the  symmetry  plane  and  directed  into  the  soil  from  the  day 
surface  [3]. 
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The  thermal  conductivity  along  the  axis  z  in  the  soil  is  negligibly  small. 

The  coordinates  of  the  pipe  center  {0,  H),  where  //  is  a  depth  of  the  pipe  laying,  and  the 

equation  of  the  pipe  circle  has  a  form:  x  +[y-H  )  ■ 

An  area  of  the  equation  integration  (2): 

|o<r<r^;  -oo<x<oo;  0<y<co;x^  +(y-H)^  >rQj  , 
with  the  bound  conditions  on  the  day  surface: 

—  cr^ir  ^  ~  ^air)  ~  ®  T  ~  ®  >  (3) 

at  the  bound  of  the  pipe  and  soil  the  conjugate  conditions  are  written,  that  are  equalities  of  the 
thermal  flows  and  temperatures  in  the  form: 

9s  “  9liq  5  ~  ^  ^liq  ’  W 

where  Ag  is  the  coefficient  of  the  soil  thermal  conductivity;  is  the  coefficient  of  a  heat 
transfer  from  the  soil  to  air;  rq  -  is  a  radius  of  a  thin-wall  pipe.  Because  the  area,  where  the 
solution  is  determined,  is  symmetrical  relatively  to  the  axis  Oy  ,  the  condition 


Ag  —  =  0  at  0  <  y  <  H  -  tq;  H  +  fq  <  y  <  CO  (5) 

dx 

is  added  at  the  bound  v  =  0  . 

The  endless  integration  area  in  the  soil  is  transformed  into  the  circle  one  using  a  conform 
transformation  [3]:  ( //  =  In  7? ) 
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H  -rq  ,  where  (»  =  exp  (//  + ;^) 


The  equation  (2)  in  new  dimensionless  coordinates  has  a  form: 
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Boundary  conditions  (3)  -  (5)  in  variables  /v,^  in  the  surface  o)  take  the  form: 
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The  symmetry  conditions  relatively  to  the  axis  Oy  : 
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Symbols  f^  ,(p  are  the  polar  coefficients  in  a  complex  plane  co  ; 
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To  construct  a  numerical  algorithms  it  is  necessary  to  express  the  derivatives  on  normal 
both  at  the  day  surface  boundary  and  at  the  pipeline  plane  by  new  variables. 

In  the  result,  derivatives  on  variables  x  ,  y  will  be  written  in  the  form: 
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On  the  day  plane  at  ju  =  0  : 
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^yy=Q 

The  derivative  on  normal  to  the  pipeline  surface  coincides  with  the  derivative  in  radius  of 


a  polar  coordinate  system  in  the  physical  surface  connected  with  the  pipeline  axis: 
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These  expressions  for  derivatives  are  necessary  to  write  the  conjugate  conditions  at  a  contact 
pipe  surface  with  the  soil  (4)  in  new  variables. 

Transformation  of  thermal  exchange  equation  (1)  at  the  flow  turbulent  regime  in  the 
pipeline  [1]  . 
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where  ©  =  1  -  —  V  = 

T  W  n,  P  W  a 

dimensionless  transformed  (deformated)  pipe  radius;  e  -  const  is  a  stretching  parameter  of  a 
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pipe  physical  radius  which  is  connected  with  the  Reynolds  dependence:  s  =  1/(0.005  Re^^^-lj. 
Using  of  the  physical  radius  transformation  is  necessary  to  compensate  a  large  near-wall 
profiles  gradient  of  a  velocity  and  temperature  at  the  turbulent  flow.  mean  and 

maximum  values  on  the  flow  velocity  axis.  P^  -  2W^rQ  I  a  ,  P^  ,  Pj  are  the  Peckle  and  Prandtl 
numbers  laminar  and  turbulent 

The  conjugate  conditions  (4)  are  written  in  the  form: 
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An  attempt  to  construct  an  algorithm  of  a  direct  solution  of  the  conjugate  problem  of  the 
pipeline  thermal  exchange  with  the  soil  was  made.  To  pursue  this  aim,  a  method  of  variable 
directions  or  the  splitting  method  on  spatial  variables  of  N.N.  Yanenko  [4]  was  applied  to 
integrate  the  equations  (6).  Two-dimensional  equation  is  splitted  into  two  one-dimensional 
equations: 
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Integrating  areas  sections  (//q  ^  M  -  O)  and  (O  <  ^  <  ;r)  are  separated  into  N-1  and  M-1 

accordingly  with  steps:  (aFq-'  =  ),  (A//„  =  //„+i  -  //„  )and  -Pm)- 

The  system  is  integrated  at  every  step  on  time.  Indexes  m ,  n  run  the  values  («  =  2,  3,  4,..., 
N-2,  N-\),  (in  =2,  3,  4,...,  M-2,  M-l)  accordingly.  The  difference  equations  for  system  (13) 
can  be  presented  in  a  three-diagonal  form: 

+  Bln  +  Cm  uCi 

^2  m+l  +  ■®2m  ^2  m  ^2m  ^2  m-l  ^  ^>2^  ■ 

The  expressions  for  the  coefficients  A,  B,  C,  D  are  written  by  the  known  form.  The  passing  ratio 
for  the  first  trinominal  equation  is  written  as: 

^In+l  ~^ln+l  C\n  +  ^1 «  • 


Using  the  boundary  condition  (7)  in  the  difference  form,  it  can  be  obtained  the  initial  values  for 
the  passing  coefficients: 
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Further,  all  passing  coefficients  ifj  „  ,  Tj  „  on  the  diminishing  index  values 
n  =  N -\,N -2,...,'i,2.  are  determined  according  to  the  re-current  formulas.  At  «  =  2  a  direct 
passing  expression  is  written 

Ui2=K,2Un+L,2.  (14) 


Therefore,  it  is  impossible  yet  to  perform  a  direct  passing  for  C/i„  determining.  Because  value 
Uii  is  unknown,  it  shall  be  determined  from  the  conjugate  conditions  (12),  but  previously,  the 
difference  trinominal  correlations  are  to  be  obtained  for  the  thermal  transfer  in  the  pipeline,  i.e.  the 
difference  scheme  for  the  equation  was  constructed  (11).  To  approximate  the  hyperbolic  left 
section,  a  four-point  “block”  scheme  was  used,  it  is  absolutely  stable  and  of  the  second  order  of 
accuracy.  An  operator  of  the  second  order  is  written  on  a  three-point  template  with  the  step 
Api  =  -  Pi  ,  where  index  /  passes  the  values  I  =  2 ,3  L-2  ,L-l .  With  increasing  of  the 

radius  p  ,  increasing  from  the  pipeline  axis  is  observed,  and  the  running  correlation  has  the  form: 

®/-l  =  +  Li  ■ 


From  the  symmetry  conditions  the  initial  values  of  the  running  coefficients  K2  =1,  L2  =0.  are 
determined  on  the  pipe  axis.  At  the  upper  bound  (on  the  pipe  wall)  the  correlation 
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(15) 

is  obtained  which  encloses  the  conjugate  conditions.  If  to  write  the  equality  (12)  in  the  form  of  a 
difference  scheme,  so  four  equations  relative  to  four  unknown  ones  are  obtained,  they  are  (12), 
(14),  (15).  Finally,  the  value  of  the  soil  temperature  at  the  contact  plane  will  take  the  form: 

+(1)Li2+ -!)(!- ^) 
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where 
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pipeline  and  solution  of  the  splitting  equation 
^  (13)  relatively  to  C/j  .  Further,  the  equation  SO lu- 
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The  algorithms  of  nonstationary  heat  exchange  were  previously  checked  in  independent 
calculations  both  for  the  soil  with  a  isothermal  pipe  and  for  the  turbulent  flow  regime  in  the 
pipe.  The  calculations  were  compared  with  results  of  [1,  3]  and  are  presented  in  Figs.  1  and  2  . 

Figure  1  shows  a  comparison  of  our  numerical  calculations  (round  points)  with  the 
calculation  results  of  [3]  under  set  conditions. 

Figure  2  presents  the  Nusselt  number  Nu  (Z,  Fq)  at  Re  =  10^,  Pr  =  1,  for  the  turbulent 
flow  in  the  pipe  and  for  uneven  change  of  the  heat  flow  2(F’q)  =  const  on  the  pipe  wall,  our 
computations  in  comparison  with  calculations  made  in  [1]  are  marked  by  a  dashed  line. 

The  calculations  of  the  certain  pipeline  are  presented  in  Fig.  3. 

ro  =  0.36  m,  H  m  ,  Re  =  10  ^  1.28  kcal/m-A  °  C,  Pr  =  1,  P  1, 


2.37  *  10^ 


m^l  h,  «air=12.8  kcal/m-h°C,  T  =  80  °  C,  t  =  15  °C; 


A  3.2  *10  kcal/M-h  C  ; 

In  the  Fig.  3  shows  dependence  of  the  mean  mass  temperature  liquid  from  distance 
and  times. 
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